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ABSTRACT 


Using the field emission retarding potential method true work 
functions have been measured for the following monocrystalline sub- 
strates: W(110), W(lll), W(100), Nb(100), Ni(100), Cu(100), Ir(110) 

and Ir(lll). The electron elastic and inelastic reflection coefficients 
from several of these surfaces have also been examined near zero . 

j 

primary beam energy. 


I. 


Work Function Measurements by the 
Field Emission R etarding Potential Method* 

by 

R.W. Strayer ' , W. Mackie, and L.W. Swanson 
Linfield Research Institute, McMinnville, Oregon 97128 

INTRODUCTION 

Of fundamental importance to the experimental and theoretical 
understanding of surfaces is the reliable measurement of its work func- 
tion. Motivated by the advancing technology of electron emitting devices 
and, more recently, by progress in experimental and theoretical under- 
standing of surfaces, an increasing number of measurements of clean, 
mono-:crystal face work function have been reported. Primarily due to 
the perennial problem of surface purity and partly due to inadequacies 
in the theoretical description of the various modes of electron emission 
one finds an inordinate degree of disagreement in the literature values 
of work function. However, the rapid advance in the methodology of 
fabricating ultra pure metals along with the increasing number of ways 
of cleaning and sensitively detecting minute concentrations of surface 
impurities is gradually eliminating surface contamination as a major 
factor in arriving at a consensus as to the appropriate clean mono- 
crystal face work function values. A further problem of lesser importance 

*This work was supported by a grant from the National Aeronautics and 
Space Administration, Lewis Research Center, Cleveland, Ohio. 

+ In partial fulfillment of Ph.D. requirement at Oregon State University. 



is that of maintaining a desired crystallographic orientation at the sur- 
face during the cleaning procedure; this problem can be detected by 
LEED and Eaue x-ray analysis and occurs for relatively few materials. 

It is in regards to the problem of the inadequacies of the existing 
models for various types of electron emission (e.g., thermionic, photo 
or field emis sion) to provide unambiguous values of mono-crystal face 
work function that has prompted this work. Here we are concerned not 
with contact potential difference measurement but rather with absolute 
(or true) work function. 

The primary electron emission processes whereby work function 
values of the emitter can be theoretically obtained are thermionic, photo 
electric and field emission. The assumptions, limitations and applic- 
ability of each of these methods for obtaining work function values have 
been adequately described in the literature. Suffice it to say that the 
Sommerfeld free electron model, which is the basic model utilized by 
the experimentalist for most electron emission processes, is severely 
strained in its application to a wide variety of non- free electron refrac- 
tory metal emitters. This limitation is succinctly described for thermionic, 
field and photo electric emission by Itskovich. Recent field emission 
energy distribution measurements have given dramatic experimental 
evidence as to the inability of the Sommerfeld based Fowler-Nordheim 


theory to explain the results from all crystallographic directions of 

2 , 3,4 

tungsten, molybdenum and copper emitters. Besides the funda- 


mental problem of model applicability, a host of minor complex effects 




such as the temperature dependence of the work function, variable re- 
flection coefficients and. the Schottky effect must be carefully incorporated 
into the theoretical framework of the emission process. 

A method of measuring the true work function of an electron col- 
lector surface which circumvents most, if not all, of the above mentioned 
difficulties is the field electron retarding potential (FERP) method. The 

FERP approach to work function studies, introduced many years ago by 
5 

Henderson , has been largely neglected with the exception of recent 

6 7 

studies of polycrystal surfaces by Holscher and Kleint. Yet this 

8 

method and an experimentally complex adaptation of the Shelton method 

are the only ways by which non-relative work functions of an electron 

collector surface can be measured. As will be shown in the following 

section, the success of the FERP method rests on the theoretical and 

experimentally verifiable fact that the voltage threshold for collection 

o 

of field emitted electrons occurs at the Fermi level E^ at 0 K or can 
be described by a Boltzman tail, i. e. , exp(E^ - E)/kT, at temperature 
T. The several experimental studies of the total energy distribution 
(TED) now in the literature^’ provide a firm basis for the validity 

of the preceding fact; also if the appropriate crystallographic direction 
of a refractory metal emitter is employed, the free electron based 
Fowler-Nordheim model of field emission is adequate for this applic- 
ation. 

In the following sections we shall describe the theoretical basis 
of the FERP method, its experimental application to the measurement 
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of mono- crystal face work functions, and the results obtained from 
several substrates. An interesting and useful fallout from the experi- 
mental approach described here is the ability to accurately measure 
the elastic and inelastic reflection coefficient for impinging electrons 
to near zero volts energy. 

THEORETICAL CONSIDERATIONS 
The expression for the differential field emitted current dl c be- 
tween energy e and e + de (where e = E - E^) in the case of a free 
electron model is given as follows: 

dl /d e =I e £/<3 / [(1 + e e/pd )d] (1) 

c o 

where p = kT/d is a dimensionless parameter. The value d .is given by 

d ='ReF/2(2m0 e ) 15 t(y) = 0. 976F/ 0 J t(y) (eV) (2) 

o 

where the electric field F and the emitter work function 0 g are in V/A 
and eV, respectively. The maximum emitter current L in Eq. (1) is 
given by the well known Fowler- Nordheim equation 

e 3 F 2 A u 

I = — exp [-4(2m0 3 ) v(y)/3fteF] 

8 it h0t (y) 

(3) 

= - 1 -- 5 ? X — F 2 Aexp [-0.683 0^v(y)/F] (A ) 

0t^(y) O 

where A Q is the area of the emitting surface from which the collected 

current originates. The image correction terms t(y) and v(y) are 

o is 1 

slowly varying tabulated functions of the auxiliary variable y = (e F) /0. 
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From Eq. (1), it is apparent the dl c /de turns on abruptly at the 
emitter Fermi level when p is small and decays exponentially with 
decreasing electron energy. The value of the half width A of the TED 
can be obtained from Eq. (1) so that at p = 0, A is given by 


A = 0.69d. (4) 

Since the practical value of d varies from 0. 1 to 0. 3 eV, the experi- 
mental half widths fall in the range 0. 07 to 0.2 eV. 

For the retarding potential method, as diagrammed in Fig. 1, the 
emitted electrons can be collected at a metal surface of work function 
0 only if their total energy E meets the condition 

E > 0 c + E f - V c 

where-. V is the emitter-to- collector bias potential; thus, decreasing 
V c allows all electrons down to the energy level e = 0 C - V c to be 
collected at 0°K. The condition V c = 0 c represents the current cutoff 
since electronic states above are not populated, and the total col- 
lected current I c at a specified value of e is given by 


A o 

Ic ~ d 


/. 


e/d 


d£ = I ( 1 - e 
o' 


e/d 


(5) 


By rewriting Eq. (5) in the working form 


lo Sl0 (I o - V"., = V 2 - 3d - V c /2. 3d (6) 
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it is clear that the values of 0 C and d can be obtained from the inter- 
cept and slope respectively of a plot of log^Q(I Q - I c )/I c versus V c . 

At emitter temperatures above 0°K log Al/I versus V deviates, 
from linearity due to the Boltzman distribution of electrons in states 
above E^. The theoretically expected effect of temperature on the TED 
has been verified experimentally and is of little consequence to the 
accuracy of utilizing Eq. (6) to obtain 0^ at T ^_300°K. This can be 
verified by noting that the temperature accounts for only a small devia- 
tion from linearity in the log Al/I 0 versus V c plot near V c = 0 C as shown 
in Fig. 2. 

Alternatively, one may obtain 0 C by noting that the value of V c 
at I c / I Q = 0. 5 when inserted into Eq. (6) yields 

0 = VJ1/2) - d In 2 (7) 

c c 


where V c (l/2) is the value of V c for which I c /I 0 = 0. 5. Eq. (7) is 
strictly applicable only for T = 0°K; however, the temperature correc- 
tion to V c (l/2) is minor and only amounts to 'v 10 meV at 300°K. The 
principal source of experimental error in this method stems from the 
uncertainty in the experimental values of d and I q due to electron 
reflection. This is discussed in greater detail in the next section. 

Eq. (1) may be differentiated with respect to e in order to obtain 


the difference in energy between the peak of the TED and the Fermi 
energy level: 


= kT In -ilL. . 
E p d-kT 




This equation, plotted in Fig. 3 at several values of d, may be used to 

obtain the theoretical value of e_ which is equal to 0_ - V . Since V_ 

P c P P 

(the position of maximum dI c /dV c on the energy axis) can be obtained 
experimentally, the value of 0 c can be obtained directly from the TED 
curve and Fig. 3. Since in practice = 30 mV this method gives 
0 c easily within 1% accuracy. 

Even though the assumptions of the Sommerfeld free electron 

model, upon which Eq. (1) and the subsequent equations are based has 

recently been found to be inadequate for certain crystallographic 

2 3 

directions of tungsten and molybdenum , the occurrence of the emis- 
sion threshold at E^ was unchanged for clean emitters. In any case, 
inadequacies in Eq. (1) due to band structure effects can be easily 
avoided for this application by choosing an emission direction (e. g. , the 
<11 1> or <310> of tungsten) for which the corresponding TED curve 
agrees well with Eq. (1). For that reason we have utilized for this 
study oriented field emitters with these directions along the emission 
axis. Thus, all that must be. known concerning the emitter in order to 
apply the FERP method of work function determination is the value of 
d which can either be calculated with sufficient accuracy from the I(V) 
characteristics of the emitter using Eq. (2) or determined experimentally 
from Eq. (6). 


ELECTRON REFLECTION 

Only one property of the collector that can detract from the com- 
plete applicability of the above equations in evaluating 0 c is electron 
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where a( e ) = 1 - R( £ ) is the electron acceptance coefficient and 
£ = | 0 c - V | is the maximum kinetic energy of the collected electrons. 
Near the collection current threshold (i.e. , e % 0) the last term of 
Eq. (9) will be small compared to the first since Ip -* 0 as e -»■ 0 and 
da /de is normally small at e = 0. 

However, as the energy of the primary electron beam increases 
above the threshold voltage considerable change may occur in a ( e ) 
(i.e., d a /de becomes large) which in turn will cause serious devia- 
tion in the apparent value of Ip. Thus, plotting the data according to 
Eq. (.6) in order to obtain an accurate value of 0 C and d will not be 
possible. In like manner it will be difficult to utilize Eq. (7) in order 
to obtain 0 C due to the inability to obtain an accurate value of I Q . 

In contrast, Eq. (8) is basically unaffected by reflection since 
the last term of Eq. (9) can usually be neglected at e = £p. We should 
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also point out that a cursory examination of Fig. 3 reveals that is 
very small (less than 40 mV) at practical values of d and T so that 
uncertainties in the exact position of £p due to reflection will not intro- 
duce appreciable error in the value of 0^. Thus, in the event that 
detectable reflection should occur for a particular collector at threshold 
the evaluation of 0 C should be accomplished from the TED curve through 
Eq. (8). 

EXPERIMENTAL APPROACH 

The basic requirement of the electron optical system for this 
application is to transform the highly divergent electron beam into a 
co- linear beam normal to the collector substrate surface and to simul- 
taneously decelerate it to zero volts. In order to maximize the analyzer 

12 

energy resolution the electron source must be highly apertured which 

in turn causes a very low beam transmission coefficient of the order of 

10" ^ to 10~^. However, if the emitter is to be operated at room temp- 

13 

erature the resolution of the analyzer need only be ^ 100 mV ; therefore, 
an electron optical system which sacrified unnecessary resolution was 
designed for this application in order to obtain a larger collector current 
to speed the data acquisition. Rather than aperture the primary beam 
to the usual ^1° half angle 0 in order to maximize resolutioiijWe chose 
for this application 0 - 8°; depending on the orientation of the emitter 
this aperture angle allowed a beam transmission of the order of 10%. 


Currents in the 10 A range were easily obtained in the focused spot 
thereby allowing the gun to be used as an electron source for other applic- 
ations as well. 

The electrostatic focusing system used in the analyzer shown in 
Fig. 4 consists of an anode, two Einzel lenses and a 500 line/in deceler- 
ating mesh electrode which established parallel equipotential surfaces in 
front of the collector. All electrodes were made from molybdenum. 

A two stage electrostatic focusing system with a virtual cross-over in 
front of the first lens was chosen over a single stage because of its 
greater optical efficiency. The lens circuitry is shown in Fig. 5, and 
the lens operating voltages given in Table I were arrived at by analytical 
computer analysis and confirmed experimentally. The anode electrode 
controls the emission level; varying the voltages on the downstream 
focusing electrodes have negligible effect upon the emission current. 

As the beam enters the first einzel lens it is partially decelerated and 
forms a virtual image of the source Zmmr behind the emitter tip. 

The second Einzel lens focuses the virtual tip image into a % 0. 5 mm 
spot size at the mesh electrode Eg. Further deceleration occurs be- 
tween electrode Eg and the mesh Eg. In most cases the mesh was oper- 
ated between 5 and 10 V relative to the collector thereby providing a 
nearly field free region between the mesh and the collector at the current 
threshold. Examination of the spatial characteristics of the beam showed 
that no significant space charge expansion of the beam occured down to 
the cut-off voltage of the mesh. By varying the screen voltage to lower 
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values and measuring the transmitted current in a Faraday cage, the 
energy distribution curve. of the electron beam passing through the mesh 
was found to be in agreement with the theoretical shape. Also, from the 
position of the current threshold the mesh work function was found to 
be approximately 4.6 eV. 

The lens system was aligned and mounted securely on four longi- 
tudinal glass rods. Both the emitter and anode could be removed as a 
unit from the tubular anode holder. In this way the emitter, which was 
held in place by a Corning 1720 glass bead in a molybdenum tube, could 
be easily replaced and prealigned in the center of the 10 mil anode 
aperture prior to insertion into the anode holder. By positioning the 
emitter in the plane of the anode aperture, no interception of the primary 
beam occurred at the anode or subsequent elements of the first lens. 
Aperturing occurred in the nearly field free region of the second Einzel 
lens by placing two 40 mil diameter stops in the last lens tube. Thus, 
electron induced desorbed ions or neutrals from the anode was elimin- 
ated and the high positive potential saddle at the anode prevented ions 
generated beyond the first aperture from bombarding the cathode. This 
design feature greatly improved the current stability without requiring 
rigorous outgassing of the electrodes. The angular convergence of the 
beam at the collector was fixed by geometry to be <1.4° for a well 
focused spot. Hence, negligible loss in resolution resulted from the 
angular deviation of the beam from perpendicularity at the collector. 
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The large aperture angle of the analyzer necessarily reduced the 

resolution of the tube as a retarding energy analyzer. Using the voltage 

separation between the 10 and 90% points on the leading edge of the 

13 

energy distribution as described by Young and Kuyatt , the resolution 
of the gun was determined to be 50 to 80 mV. This resolution was 
adequate to resolve the leading edge of the energy distribution curve 
at room temperature. 

The emitter orientations selected for the field electron source 

in this study were <1 1 1> and <310> tungsten, fabricated from zone 1 ' 

2 

oriented wire. Previous studies have shown that electrons field 
emitted from these orientations exhibit energy distribution curves 
that agree closely with the Scrrmo-feld free electron model upon which 
the theoretical expressions of the previous section are based. Further- 
more, the work functions of the crystal planes intersecting these direc- 
tions are quite low -- 4.3 eV for the (310) plane and 4.4 eV for the (111) 
plane -- thereby providing the highest beam transmission values. 

The single crystal collector substrates of this study were shaped 
and mounted in the holder as shown in Fig. 4. The face of the collector 
crystal was circular with a diameter of 200 mils. This was sufficiently 
large compared to the 20 to 40 mil beam size to eliminate edge effects. 
Thermal and electron induced desorption cleaning of impurities at the 
collector crystal -were accomplished through electron bombardment. 
Collector crystals could be easily replaced by removing the glass seal 
which holds the collector support rod. 
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EXPERIMENTAL PROCEDURE 

The single crystal collector surfaces were fabricated from MARZ 
grade material by a high speed grinding wheel and electro- chemical 
machining techniques. Several mils of the collector surface were re- 
moved by electrochemical etching in order to eliminate mechanical 
defects. The etchant solutions for the nickel, copper, tungsten, niobium 
and iridium surfaces were concentrated phosphoric acid for both nickel 
and copper, sodium hydroxide, 25% hydrofluoric /25% sulfuric acid 
solution in water, and 5% sodium hypochlorite, respectively. The 
alignment of the desired crystal directions with respect to target normal 
was within ±1° as shown by Laue x- ray examination. Monocrystallinity 
of the substrates was - carefully checked both before and after measure- 
ment by a high powered optical microscope and Laue x-ray examination. 

After mounting the collector substrate, the tube was evacuated 
to vlO ^ torr pressure range. The arrangement allowed for the option 
of immersing the tube partially or completely in liquid nitrogen in order 
to enhance vacuum stability and to extend the temperature range of the 
work function measurements. 

Cleaning of the crystal surface was accomplished by electron 

o o 

bombardment heating to 2100 K for niobium and tungsten, 1700 K in 

the case of iridium, 1400°K for nickel and 1100°Kfor copper. The 

copper crystal was exposed to 10" ^ torr of hydrogen and heated to 



1250 K. In order to remove potential carbon contamination the nickel 

-5 

crystal was first heated in 10 torr of oxygen followed by heating in 
10 torr of hydrogen and then heated to ^ 1400°K in high vacuum. 
Thermal heating was continued until the field electron emitter, which 
was also cleaned thermally, and therefore very sensitive to gas release 
from the collector, showed no change in current after flashing the col- 
lector substrate to its cleaning temperature. As a corollary check, 
the absence of further change in the collector work function on heating 
was used as an indication of a clean surface. 

The emitter-to-collector current- voltage characteristics were 
taken several times for each collector crystal and plotted on an x-y 
recorder. A computer program was formulated to plot the data accord- 
ing to Eq. (6) so that a value of 0 c and d could be obtained. As will 
be noted later in the paper this method of determining 0 c and d was 
not always applicable due to reflection. Therefore, the differential 
curve was also taken by utilizing the position 3 circuitry shown in Fig. 5, 
which involves the well known electronic differentiation method utilizing 
a PAR HR- 8 Lock-In Amplifier. A 10 mV 1000 Hz signal v s was fed 
to the emitter through the transformer T-^. The magnitude of this sig- 
nal detected by the Lock-In Amplifier (LA) is proportional to dI c /dV. 
From the position of the TED peak and Eq. (8), a-value of 0 C was ob- 
tained; this result could be obtained with an experimental accuracy of 
± 20 mV. The I- V and TED data were usually taken at emitter temp- 
eratures of 77 and 300 °. K. 



EXPERIMENTAL RESULTS 


The values of 0 C obtained by the FERP method are reported in 
Table II along with a comparison with other methods. These values 
were obtained via Eq. (8) from the TED curves shown in Figs. 6- 8. An 
interesting feature of these results is the additional structure in the 
TED for Cu(100), Nb(100) and to a lesser degree for Ir (111). This 
anomalous structure near the threshold of the TED curve was not as 
apparent in the other results and was determined to be due to unusual 
electron reflection. The integral current- voltage characteristics shown 
in Figs. 9-16 on a compressed voltage scale clearly show that the 
origin of the TED structure in Cu(100) and Ir ( 111) is due to electron 
reflection which varies rapidly at the threshold voltage. Although the 
integral curve for Nb(100) was not obtained due to inadvertent melting 
of the crystal, the structure in the TED curve suggest a reflection co- 
efficient which also varies rapidly near threshold. All other substrates 
show reflection at threshold, but to a lesser degree. 

If we assume that all the reflected electrons are collected at the 

mesh the sum of the collector current I and mesh current I is given 

c s 

by 


I' 

P 


I + I 


(10) 


where 1^ is the emitter current arriving at the mesh. Noting that the 
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mesh transmission is given by I^/I^ = T, where 1^ is the current imping- 
ing on the collector, and that (1-R) = I c /Ip, one obtains 


(1-R) = I c /T(I c + I a ). 


(ID 


Since T, I_ and I are measurable quantities R may be determined as 

a function of V c as shown in Figs. 9- 16 for each of the substrates. It 

is important to point out that Eq. (11) is valid when 1^ is at its saturated 

value. Since we do not wish to count electrons reflected by the retarding 

field, it is necessary to use the theoretical variation of I with V to 

1 pc 

calculate R in the retarding range. As a self consistent check on the 
electron current accounting we can also measure 1^ in the deep retard- 
ing range preceding the threshold where I ' = I g ; this value of I ' generally 
agreed with that calculated from Eq. (10) at- saturation. The theoretical 
dependence of the current I on V is also plotted in Figs. 9-16 with 

p c 

appropriate normalization factors applied. Clearly the values of R for 
both Ir(lll) and Cu(100) are unusually high and rapidly changing near 
threshold. It must be remembered that the energy spread of the primary 
beam at half height is % 0. 2 eV so that structure in R less than 0. 2 eV 
width may be detected but not accurately reproduced with respect to 
shape. However, with these definitions in mind, one may confidently 
measure R down to zero primary beam energy. It should be pointed 
out that in order to obtain the value of primary beam energy 0^ must 
be subtracted from the abscissa of each of the I-V curves. The latter 
correction has been made in Figs. 9-16. 
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In our experimental setup when V c - 0^ > V s - (s refers to 
the screen-mesh electrode) those reflected electrons which lose energy- 
through inelastic processes will not be collected at the mesh but instead 
will be returned to the collector. Thus, by fixing V g = 5.8 V only 
specularly reflected elastically scattered electrons will escape from the 
collector crystal. The elastically reflected electron coefficient R g has 
been measured for several substrates as given in Figs. 10-15. By setting 
V = 130 V all reflected electrons in the energy range investigated return 
to the mesh and total reflection R t curves are obtained (see Fig. 10). 
Hence, it is possible by this technique to measure the inelastically re- 
flected electron coefficient R^ n by noting that R- n = R^. - R g . In Figs. 

11-15 values of R. are given for several of the crystal faces. For the 
in 

crystal faces ;Cu(100) and W(lll) only R fc values were obtained. Diffi- 
culty with leakage current to the mesh electrode caused some minor 
uncertainty regarding the absolute values of R for Ir(110), W(110) and 
Cu(100). 

DISCUSSION 

Work Function 

As shown in Fig. 6 and Table II the work function values of 4.47, 

4. 63 and 5.25 eV obtained for the (111), (100) and (110) planes of tung- 
sten are in close agreement with values reported by other methods, 
particularly thermionic and field electron emission methods. Since 
the electron reflection coefficient is low (less than 15%) near threshold 
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for the tungsten results, the values of 0 C obtained from Eq. (6) compare 
within experimental error with the values computed from the TED curves 
via Eq. (8) (see Fig. 2 for a typical result utilizing Eq. (6)). This 
merely provides the expected self consistent check on the experimental 
method when R is small or a slowly varying function of energy near 
threshold. 

Notable, is the well known discrepancy between 0(110) obtained on 

macroscopic crystals and microscopic crystals employed in field emission 

(i. e. , by use of the Fowler- Nor dheim (FN) equation) techniques. This 

discrepancy is believed to be due to large geometric facets which occur 

on W(110) planes of a field emitter and which can alter the FN results 

in the observed fashion. Tungsten work function values exhibit the usual 

increase with surface atom density in accordance with the semiquantitative 

2 1 

theoretical expectations of: Smoluchowski. 

The Ir(lll) plane, possessing one of the highest atom densities, 

yields the expected high work function value 0 = 5. 76 eV as shown in 

Fig. 8. A similar high value of 0 for Ir(lll),: in remarkable agreement 

15 

with the value obtained here, has been reported by Zandberg who ob- 
serves this to be the highest work function value reported for a macro- 
scopic single crystal metal surface. The agreement is surprising in 
view of the initial polycrystalline nature of Zandberg 's iridium crystal. 
Our results indicate that this crystal face is relatively stable with respect 

to thermal induced recrystallization and faceting effects such as observed 

19 

on the high work function Re(0001) face. 
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The Ir(llO) work function is somewhat lower than Ir (111) as expected 

»v- ... 

on the basis of the lower atom density. Considering that the (110) face 
of Ir is the third most densely pack plane, the relatively large value 
0 = 5.42 eV is indicative of the generally high work function of all crystal 
faces of Ir. 

The value of 0 = 4. 18 eV for the (100) face of bcc niobium shown in 

16, 

Fig. 7 agrees closely with thermionic and FN values reported by others. 

(See Table II). The work function values for W(100) and Mo(100) obtained 

by a variety of methods indicate best values of 4.65 and 4.40 eV respec- 
17 

tively. Both W(100) and Mo(100) possess identical atom densities 
14 2 

cr = 10.9 x 10 atoms/cm . Since Mo and Nb are neighboring elements 

in the periodic table one might expect geometric factors to dominate there- 

' > 

by causing a larger 0 for Nb. This is obviously not the case. These 

results underscore the role of electronic factors in the variation of work 

function with crystal face and material. . 

Few reliable values of work function for the (100) face of Ni and Cu 

are reported in the literature - obviously, due to the surface cleaning 

difficulties. Our results, shown in Fig. 7 and Table II are in rough agree- 

18 20 

ment with photo electric and thermionic values reported for Ni(100) 

and Cu(100) respectively. Using careful cleaning procedures and a gold 

22 

reference electrode, Delchar reports a work function value of 5. 16 eV 
for Cu(100) in excellent agreement with this work. Delchar points out 
that any possible error in the work function value for the reference elec- 
trode would likely be in such a direction as to raise the Cu(100) value. 


Our value of Ni(100) work function is 0.30 eV larger than other 
reliable measurements.- Most previous measurements on crystals of 
questionable surface purity converge on a value of 4.9 eV. Clearly, a 
combination of the FERP method with Auger analysis would be extremely 
helpful in answering the question of surface cleanliness. 

It is interesting to note that the elements Cu and Mo are each 
adjacent to Ni and Nb respectively in the periodic table. In each case 
the change in the free atom electronic structure between adjacent elements 
involves the filling the inner d orbitals. In the case of the fee metals 
Ni and Cu the work function of the (100) crystal face decreases with 
atomic number Z, whereas with the bcc metals Nb and Mo the work 
function of the (100) face increases with Z. At the same time each of 
these pairs of metals exhibit a slight increase in their respective lattice 
constants with Z. The value of the (100) work function for the bcc metal 
W is larger than the Mo value by M). 25 eV even though both metals 
possess nearly identical lattice constants and differ in electronic struc- 
ture by the filling of inner 5p and 4f orbitals. Clearly, a predictive 
theory of metallic work function must include complex structural and 
electronic factors. 

Electron reflection 

■ The simultaneous measurement of work function and electron re- 
flection was undertaken here primarily as a further index of surface 
cleanliness. The ease of separating the elastic R g and inelastic R^ n 
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reflection coefficients by this method motivated a more detailed study of 
these coefficients as a function of the primary electron energy Ep. Several 
interesting and surprising features have been observed in the variation of 
R with Ep in the low energy range. 

Fig. 9 shows that the variation of R^ with Ep for W(lll) is in good 

23 

agreement with careful measurements of R g by Armstrong. Because 
of the ability of the FERP method to accurately measure R to within a 
few tenths eV of threshold our results show a definite peak at 3 eV here- 
tofore unobserved. A small peak observed at 6 eV is in agreement with 
Armstrong's results. Disagreement in absolute scale at larger values of 
Ep can be ascribed to the difference between Armstrong's measurement 
of R e and in this work the measurement of R t for this crystal face. As 
also shown in Armstrong's results surface contamination markedly alters 
R 0 ; thus the close agreement between the two results mutually supports 
the claim of an atomically clean surface. 

The W(100) results in Figs. 10 and 11 show the experimental vari- 
ation of I with Ep for V <3 = 100 and 5.25 V respectively, thereby illu- 
strating the difference between the total and elastic reflection modes. In 
Fig. 11 the value of the experimental reflection coefficients R- n and R g 
are given along with Armstrong's results for R g . Both results agree 
that a sharp peak in R g occurs between 3. 5 and 4. 0 eV as observed simi- 
larly in the W(lll) resultSi In disagreement with Armstrong's results, 
the large peak in R g at 8. 0 eV is missing and instead a smaller peak at 
6.5 eV appears. In addition, a major peak observed in our results, but 
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absent in Armstrong's , occurs at 16.5 eV . Both investigations agree on 
a broad peak in R e at 10. 5 eV. Differences in the R e vs Ep curves obtained 
in the two investigations must be ascribed to either a small degree of sur- 
face contamination,crystal inhomogenity or experimental artifact. 

The variation of R e with Ep for W(110) shown in Fig. 12 over a wider 
Ep range exhibits a large peak modulated with fine structure in the 4-6 eV 
range. Lesser peaks occur at 14 and 27 eV as R g decreases from a peak 

value of 0. 34 at 5 eV to less than 0.03 at 75 eV. The W(110) R g curve 

24 

agrees closely with an earlier result by Kahn, et. al. 

The inelastic reflection coefficients R^ n given in Figs. 11 and 12 for 

W(100) and "W(110) exhibit very similar structure for both crystal faces. 

That is, each exhibits an onset threshold at 3 eV, fine structure in the 

3 to 10 eV range, and additional peaks at 12. and 17 eV. The larger range 

of E for W(110) indicates an increase in R . to a value of 0.47 at 52 eV 

p i n 

where it then begins a slight drop. The R 0 and R^ n curves for W(100) 

and W(110) both cross in the 12 to 15 eV range and with increasing Ep the 

ratio R. /R increases to a value of ^ 90 at E„ = 60 eV for W(110). 
in e p 

One of the most unusual results of this study is the exceedingly 

large value of R = 0.75 as E -* 0 as shown in Fig. 13 for Ir( 1 1 1). The 
e p 

Ir(110) results in Fig. 14 measured over a larger Ep range also show an 
increase in R g at threshold but not nearly as large a value as for Ir (111). 
From a dynamical point of -view, these results suggest the occurrence of 
a band gap relatively free of surface states approximately 5.79 eV above 

24 

^Earlier results by Kahn, Hobson and Armstrong show a definite peak 
in this energy range. 



the Fermi level along the <11 lv> direction of Ir. Both crystals shpw one 
additional large peak in R g at 11.5 and 15 eV for Ir (110) and Ir (111) re- 
spectively. While, the variation of R with E differs sharply for the two 

e p 

crystal faces, the, variation of R_ with E^ is quite similar. On both cry- 
stal faces R^ n increases monotonically with E showing a series of small 
peaks in the 8 to 30 eV range. No maximum in R^ n is apparent in this 
range of E^ as was the case for W(110) and Ni(100). In contrast with the 
W results the threshold for R^ n occurs at a higher voltage, i.e., in the 
6 to 8 eV range. It is interesting to note that whereas R^ is large for 
Ir(lll) at Ep = 0, the value of R fc for Ir(110) is smaller at Ep = 0 and 
increases to a value of 0.7 at 70 eV where is appears to be still increase 
ing with Ep. 

The variation of R g with E^ for Ni(100) is in general agreement with 

diffraction intensity studies. The peaks in R g at threshold 15, 28 and 38 eV 
, ■ 25 

have been observed elsewhere, but with somewhat differing relative mag- 
nitudes. The peak near threshold can be ascribed to the first order pri- 
mary Bragg peak of the 00 band. We also observe an increase in the 
elastic peak widths and decrease in the elastic peak heights as R^ n increases 

in accordance with certain theoretical viewpoints of the electron scattering 
25,26,27 

mechanism. A large peak normally observed in R g at 55 eV is 

missing in the present results. Part of the discrepancy may be due 
to the presence of a few - grain boundaries which developed near the 
edge of our Ni(100) crystal prior to obtaining the reflection 
curves. Although the central portion of the crystal 
(where the electron beam impinged) appeared to be monocrystalline 



(as determined by Laue x-ray), a slight reduction in 0 occurred indicating 
the possibility of structural effects on the reflection results. 

The inelastic reflection threshold for Ni(100) occurs a few volts above 
zero and exceeds R g for Ep > 16 eV. A small peak in Reoccurs at 18 eV 
followed by a major maximum at 40 eV. It is clear from the Ni(100) 
results along with the above mentioned results that inelastic processes be- 
come important above 15 to 20 eV and must therefore be considered care- 
fully in theoretical attempts to explain LEED intensity variations with E^. 
The Cu(100) results shown in Fig. 16 were obtained for R^ only over 

a limited range of E . Since R. is likely to be small in this energy range 

pm 

R fc =s Rg. Although not as large as the threshold peak for Ir(lll), the 

28 

threshold peak for Cu(100) is quite large and has been attributed to the 

29 

well known band gap in the bulk E(k) diagram occuring at the vacuum 
level. Previous LEED studies of the specular reflectivity I from Cu(100) 
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gave a value of R g = 0. 36 around Ep = 1 eV. ““ The results reported here 
agree closely with the latter results but also show a larger peak at E = 0 

r 

of R = 0.47. In addition, a very small peak is observed at E„ = 2. 0 eV. 
e ■ r 

These extremely narrow peaks (less than 0. 5 eV) point out the resolving 
power of this technique near threshold and agree with the previous observa- 
tion that peak widths in R e below the onset of inelastic reflection are narrower. 
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SUMMARY 


The use of the FERP method to accurately measure true work func- 
tion of a variety of single crystal surfaces (see Table II) with a minimum 
of assumptions regarding the nature of the electron emission process has 
been demonstrated. A combination of the FERP and LEED /Auger should 
provide a powerful combination of techniques to accurately and conveniently 
measure bare surface true work functions of conducting surfaces. The 
occurrence of electron reflection from the collecting surface can generally 
be tolerated without introducing appreciable error (ca ± 0. 02 eV) pro- 

vided that the work function is obtained from the peak in the TED curve. 

An unusually large work function of 5. 76 eV for Ir(lll) was obtained in 

15 

agreement with an earlier result. Other results agreed favorably with 

previously measured and accepted values of work function. 

23 

As shown by Armstrong surface contamination or crystal imper- 
fection can alter the elastic reflectivity vs. primary energy curves. 

We have shown that the FERP technique can yield both elastic and inelastic 
reflection curves accurately to within a few tenths eV of zero primary 
energy and thereby can be used as an indicator of surface conditions. 
Values of R g and R^ n for several crystal faces agreed reasonably well with 

Of 

published curves. We noted that R^ n >> R g as E > 15 eV for most 

P 

crystal faces examined. An unusually large value of R 0 was obtained 
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for Ir (111) and Cu (100) at threshold. Whereas the structure in the R e 

curves wase highly dependent on crystal orientation, the R curves were 

in 

generally -unaltered by crystal orientation including the threshold values 

of E which were usually 3 to 8 eV. 

P 

Finally, we anticipate that the FERP technique can be profitably 

used to examine the combined effect of chemisorption on 0, R and R. . 

r e m . 
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Table L» Operating voltages on lens elements 
(see Fig. 5 with emitter at 0 V). 


Lens element 


E 5 = E 7 
E 6 = E 8 


Operating voltage 

0. 06E^ 

0 

0. 04E 3 

5-10 or 130-150 V 


27 



Table H 


Work Function (eV) 


Material 

Atom Density 
(atoms/ cm ) 

(this work) 

Thermionic Photoelectric 

Field 

Emission 

(emitter 

values) 

W(100) 

14. 1 

5.25 ± 0. 02 

14 

5.35 ±0.05 

5.9 ± 0. I 2 

W(100) 

10.0 

4.63 ± 0. 02 

. , A 14 

4.60 ± 0.05 

4.70 ± 0. 05 2 

W(lll) 

5.77 

4.47 ± 0. 02 

4.40 ± 0. 02 14 

4.45 ± 0. 03 2 

Ir ( 111) 

15. 8 

5.76 ± 0. 04 

5. 79 ± 0. 03 15 


Ir'(llO) 

9.7 

5.42 ± 0. 02 



Nb(100) 

10.9 

4. 18 ± 0. 02 

3. 95 ± 0.03 16 

3.87 ± 0. 01 17 

Ni(lOO) 

16. 15 

5. 53 ± 0.05 

18 

5. 22 ±0.04 


Cv.flOO) 

15.4 

5. 10 ± 0. 05 

4.9 20 
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, FIGURE CAPTIONS 

Fig. 1. Potential energy diagram for a field electron retarding potential 
analyzer. When the collector is biased such that only electrons 
from the Fermi level of the emitter can reach the collector, the 
battery voltage V c is equal to the collector work function 0 c> 

Fig. 2. Plot of In AI/I 0 vs. collector voltage for a W(110) collector. The 

intercept of the abscissa at AI/X 0 = 1 gives 0 C according to Eq. (6). 

Fig. 3. The difference £p in energy between the peak of the TED curve 
and the Fermi level as a function of T and energy parameter d. 

Fig. 4. Diagram of FERP tube showing pertinent features of the electron 
optical system and collector single crystal. 

Fig. 5. Diagram of the electrical circuitry associated with the FERP tube. 

Fig. 6. TED curves for W(lll), W(100) and W(110). 

Fig. 7. -TED curves for Cu(100), Ni(100) and Nb(100). 

Fig. 8. TED curves for Ir(lll) and Ir(110). 

Fig. 9. -Experimental and theoretical integral I (V c ) curves obtained from 
W(lll). Solid line shows the total reflection curve R^.. Dashed 
curve gives the experimental elastic reflection curve R g obtained 
by Armstrong. ^ 
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Fig. 10. Experimental and theoretical I (V ) curves obtained for W(100). 

Solid line shows the total reflection R^. curve. The mesh voltage 
V s = 100 V. 

Fig. 11. Experimental and theoretical I (V c ) curves obtained from W(100). 

Solid line curves show inelastic R. and elastic R„ reflection 

m e 

coefficients. The mesh voltage Vg = 5.25 V. 


Fig. 12. Experimental and theoretical I (V ) curves obtained from W(110). 
Solid line curves show inelastic R^ n and elastic R e reflection 
coefficients. <■ 


Fig. 13. Experimental and theoretical I (V ) curves obtained from Ir(lll). 

Solid line curves show inelastic R. and elastic R reflection 

m e 

coefficients. 


d 

Fi'gc, 14. Experimental and theoretical I (V ) curves obtained from If (100). 

Solid line curves show inelastic R. and elastic R reflection 

in e 

coefficients. 


Fig. 15. Experimental and theoretical I((V c ) curves obtained from Ni(100). 

Solid line curves show inelastic R- and elastic R„ reflection 

in e 

coefficients. 


d 

Fig. 16. Experimental and theoretical I (V ) curves obtained from Cu(100). 
Solid line curve shows the total R^. reflection coefficient. 



Fig. 2 



34 








Fig. 5 
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